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Abstract 

Charged particle production has been studied in neutral current deep inelas- 
tic ep scattering with the ZEUS detector at HERA using an integrated lumi- 
nosity of 0.44 fb" 1 . Distributions of scaled momenta in the Breit frame are 
presented for particles in the current fragmentation region. The evolution of 
these spectra with the photon virtuality, Q 2 , is described in the kinematic re- 
gion 10 < Q 2 < 41000 GeV 2 . Next-to-leading-order and modified leading-log- 
approximation QCD calculations as well as predictions from Monte Carlo models 
are compared to the data. The results are also compared to e + e~ annihilation 
data. The dependences of the pseudorapidity distribution of the particles on Q 2 
and on the energy in the 7p system, W, are presented and interpreted in the 
context of the hypothesis of limiting fragmentation. 
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1 Introduction 



Quark fragmentation has previously been studied experimentally in deep inelastic ep 
scattering (DIS) at HERA using observables such as multiplicity moments, scaled mo- 
mentum distributions and fragmentation functions [TH3]. The results were compared to 
those obtained in e + e~ and pp collisions. In general, universal behaviour has been es- 
tablished and scaling violations of the fragmentation functions jl] observed. It has also 
been observed that perturbative Quantum Chromodynamics (pQCD) calculations using 
the modified leading-log-approximation (MLLA) [5JE] and assuming local parton-hadron 
duality (LPHD) [7J do not provide a full description of the data. 

In this paper, multiplicity distributions of charged hadrons in the current region in the 
Breit framed are presented as functions of the virtuality of the exchanged boson, Q 2 per 
unit of the scaled momentum, x p = 2PB rc[t /Q, and the variable ln(l/x p ) in bins of Q 2 . 
Here, Pereit denotes the momentum of a hadron in the Breit frame. The data sample 
collected with the ZEUS detector between 1996-2007, comprising 0.44 fb _1 , enables the 
study to be extended to Q 2 as high as 41 000 GeV 2 . Predictions from next-to-leading-order 
(NLO) QCD calculations that combine full NLO matrix elements with fragmentation 
functions (FF) obtained from fits to e + e~ annihilation data [5HTT] are compared to the 
measurements. Predictions from MLLA+LPHD [6l[T2j[T3] and leading-order plus parton- 
shower Monte Carlo programs are also considered. 

In addition, the measurements are compared to previous ep results [IH3] and to e + e~ 
annihilation data [2]. The hadronisation in the current region in the Breit frame in ep 
scattering can be compared directly to the hadronisation in one hemisphere of e + e~ anni- 
hilation events. There, particle momenta are scaled to half of the centre-of-mass energy, 
E* = y/s/2. Previous studies on DIS hadronisation [TH5jlTo] have shown good agreement 
with e + e~ annihilation at medium Q 2 . At lower Q 2 , Q 2 < 40 GeV 2 , the two processes were 
observed to behave differently [2]. This can be explained by higher-order QCD processes 
such as boson-gluon fusion (BGF) and initial-state QCD Compton radiation occurring as 
part of the hard interaction in ep scattering but not in e + e~ annihilation. 

Finally, the density of charged particles is studied as a function of the particle pseudora- 
pidity, 77 Breit ; in bins of Q 2 , and as a function of the total 7*p centre-of-mass energy, W. 
The data are used to test the hypothesis of limiting fragmentation [16], in which there 
has recently been renewed interest, most notably in relativistic heavy ion collisions |17j . 



The Breit frame is defined as the frame in which the four-vector of the exchanged photon becomes 
(0,0,0,-Q). 
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2 Experimental setup 



A detailed description of the ZEUS detector can be found elsewhere [18]. A brief out- 
line of the components most relevant for this analysis is given below. Charged particles 
were tracked in the central tracking detector (CTD) |19], which operated in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. The CTD consisted of 72 
cylindrical drift-chamber layers, organized in nine superlayers covering the polar-angle^! 
region 15° < 9 < 164° and the radial range from 18.2 cm to 79.4 cm. Before the 2003- 
2007 running period, the ZEUS tracking system was upgraded with a silicon microvertex 
detector (MVD) [20]. 

The high- resolution uranium-scintillator calorimeter (CAL) [2T] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers 
and longitudinally into one electromagnetic section (EMC) and either one (in RCAL) 
or two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the 
CAL was called a cell. Under test-beam conditions, the CAL single-particle relative energy 
resolutions were (Te/E = 0.18 /\f~E for electrons and &e/E = 0.35 /y/~E for hadrons, with 
E in GeV. 

The position of the scattered electron was determined by combining information from the 
CAL, the small-angle rear tracking detector (SRTD) [22] and the presampler (PRES) |23j . 
both mounted on the face of the RCAL. 

The luminosity was measured using the Bethe-Heitler process ep — > eyp by a luminosity 
detector which consisted of a lead-scintillator [241 - 126] calorimeter and, after 2002, an 
independent magnetic spectrometer [27]. The fractional systematic uncertainty on the 
measured luminosity was 2% and 2.6% for the 1996-2000 and 2004-2007 running periods, 
respectively. 



3 Event sample 

The data presented here were collected with the ZEUS detector at HERA between 1996 
and 2007 and correspond to an integrated luminosity of 0.44 fb" 1 . During 1995-97 (1998- 
2007) HERA operated with protons with an energy of E p = 820 GeV (920 GeV) and 



2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 



2 



electrons^ with an energy of E e = 27.5 GeV, resulting in a centre-of-mass energy of 
^ = 300 GeV (318 GeV). 

A three-level trigger system [T51I2S1I2H] was used to select events online. It relied on the 
presence of an energy deposition in the CAL compatible with that of a scattered electron. 
At the third level, an electron, identified using the pattern of energy deposits in the 
CAL [30] and having an energy larger than 4 GeV, was required. 

The kinematic variables, Q 2 and the Bjorken scaling variable, x, as well as the boost 
vector to the Breit frame were reconstructed using the double angle (DA) method [31] 
based on the angles of the scattered electron and of the hadronic system. The total energy 
of the 7p-system, W, was calculated using W 2 = Qda(1 — ^da)/^da- 

The tracks used in the analysis had to be associated with the primary interaction vertex 
and were required to be in the region of high CTD acceptance, \r)\ < 1.75, where t] = 
— ln(tan#/2) is the pseudorapidity of the track in the laboratory frame with 9 being the 
polar angle of the measured track with respect to the proton direction. The tracks had 
to pass through at least three CTD superlayers and were required to have a transverse 
momentum, P^ rack > 150 MeV. The details of the event reconstruction are similar to those 
in a previous ZEUS publication [2] and are described in detail elsewhere [52"] . 

To reconstruct x p , the momentum four- vector of each track was boosted to the Breit frame 
assuming the particle to have the pion mass. 

The analysis of the scaled momenta was restricted to events with Q 2 > 160 GeV 2 . A well 
reconstructed neutral current DIS sample was selected by requiring the following: 

• E' e > 10 GeV, where E' e is the energy of the scattered electron; 

• He < 0.95, where y e is the inelasticity, y = Q 2 /sx, estimated from the energy and angle 
of the scattered electron; this reduces the photoproduction background; 

• 2/jb > 0.04, where yjB is estimated by the Jacquet-Blondel method [33]; this rejects 
events for which the DA method gives a poor reconstruction; 

• 35 < <5 < 60 GeV, where 5 = XX — Pz t ) and Ei is the energy of the i-th calorimeter 
cell, Pz i is its momentum along the Z axis and the sum runs over all cells; this removes 
photoproduction and events with initial-state radiation; 

• |Z vcrtcx | < 50 cm, where Z veTtex is the Z component of the position of the primary 
interaction vertex; this reduces background from events not originating from ep colli- 
sions; 

• the position (X, Y) of the scattered electron candidates in the RCAL was required to 
satisfy y/X 2 + Y 2 > 35 cm. 

3 The term "electron" is used for both electrons and positrons. 
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The limiting-fragmentation analysis was extended to events with lower Q 2 values, Q 2 > 
10 GeV 2 , using data collected during 1996-2000 corresponding to an integrated luminos- 
ity of 77 pb _1 . The same selection as described above was used but, for events with 
10 < Q 2 < 160 GeV 2 , the last requirement was modified and one additional requirement 
introduced: 

• depending on experimental conditions, the position (X, Y) of the scattered electron 
candidate in the RCAL was required to satisfy \X\ > 12 cm and \Y\ > 6 cm or 
|X| > 14 cm and \Y\ > 14 cm or \/X 2 + Y 2 > 35 cm. 

• ?7max > 3.2, where ?7 max is the pseudorapidity in the laboratory frame of the most 
forward energy deposit with at least 400 MeV; this removes diffractive events. Even 
before the cut, for Q 2 > 100 GeV 2 , the contamination with diffractive events is well 
below 5%, and the resulting corrections to the multiplicity are below 2%. In the lowest 
Q 2 bin, the contribution of diffractive events would grow up to 10% without the cut. 

Only tracks with p|* reit < enter the analysis of the scaled momenta, whereas this 
restriction does not apply to the limiting-fragmentation studies. The total number of 
events is listed in Table [TJ 

4 Models and Monte Carlo Simulations 

The NLO perturbative QCD calculations considered, combine the full NLO matrix ele- 
ments with NLO fragmentation functions obtained from fits to e + e~ data [8HTT]. The 
resulting predictions were obtained using the Cyclops program |34j . 

The MLLA calculations [6| fT2l 1351437] describe parton production in terms of a shower 
evolution. They depend on two parameters only, the effective QCD scale, A eff , and the 
infrared cutoff scale, Q , at which the parton cascade is stopped. The calculations in- 
trinsically include colour coherence and gluon interference effects. Leading collinear and 
infrared singularities are removed and energy-momentum conservation is obeyed. 

To connect predictions at the parton level to the hadron-level data, LPHD [7] was as- 
sumed, which leaves only one free parameter, the hadronisation constant, K^. The con- 
version from energy to momentum spectra for the final-state hadrons was performed 
assuming an effective hadron mass, m e fj = Q Q [13]. The input parameters for the calcu- 
lations were obtained by fits to LEP e + e~ data. Conservative uncertainties, equivalent 
to three standard deviations of the experimental uncertainty, were assumed for the pa- 
rameters: Q = A cS = 270 ± 20 MeV and K h = 1.31 ± 0.03. The A eff value agrees with 
the value A e g- = 275 ± 4(stat.)^g(syst.) MeV deduced from a ZEUS analysis of scaled 
momenta in dijet photoproduction [3B|. The usage of input parameters deduced from 
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LEP data is justified by the assumed equivalence of hadronisation in one hemisphere of 
e + e~ annihilation and in the current region of ep interactions in the Breit frame. 

The predictions from several Monte Carlo (MC) models were compared to the measure- 
ments. Neutral current DIS events were generated using the leading-order QCD Ari- 
adne 4.12 program [39]. The QCD cascade was simulated using the colour-dipole model 
(CDM) [40] inside Ariadne. Additional samples were generated with the MEPS model 
of Lepto 6.5 [41]. Both MC programs, Ariadne and Lepto, were also used to calculate 
detector acceptances and to correct the data to the hadron level. For this purpose they are 
used with the Djangoh 1.1 [42] interface and QED radiative effects are included using 
the Heracles 4.6.1 [43] program. Both MC programs use the Lund string model [44] for 
hadronisation as implemented in Jetset 7.4 [4o"jl4"6]. Hadrons are assumed stable if their 
lifetime is larger than 3 x 10~ n s; their decay products are not considered. This excludes 
in particular the charged decay products of K° and A particles. 

All generated events were passed through the ZEUS detector- and trigger-simulation pro- 
grams which are based on GEANT 3 [47]. They were reconstructed and analysed by the 
same program chain as the data. 

5 Correction procedure and systematic uncertainties 

All measured distributions were corrected to the hadron level. The correction factors 
were calculated bin by bin using MC events. For the scaled momentum spectra, the 
factors are typically below 1.2 for Q 2 < 5000 GeV 2 , but rise up to 1.5 for higher Q . 
The corrections for charged-particle densities as a function of the pseudorapidity are of a 
similar magnitude and approach 1.5 as ^ Brcit increases towards the most positive values 
measured. 

Cross sections, measured separately for each data-taking period, were combined using a 
standard weighted average [48]. The dependence of the scaled momentum distributions 
on the variation of the proton-beam energy in the different data samples was determined 
using MC events; the resulting changes were found to be smaller than 0.5% and were 
neglected. Finally, the results were corrected to the QED Born level using correction 
factors obtained from MC, reducing the charged hadron multiplicities by up to 4%. 

The systematic uncertainties were investigated separately for data with Q 2 above and 
below 160 GeV 2 . The numbers in parentheses correspond to the largest variations ob- 
served in the scaled momentum spectra. The uncertainties in the limiting-fragmentation 
distributions are of similar magnitude. For data with Q 2 > 160 GeV 2 , the systematic 
uncertainties are due to: 
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• imperfections in the simulation affecting the determination of the efficiency of event re- 
construction and event selection (tl%)- This was evaluated by modifying the selection 
cuts within the experimental resolutions. 

• an uncertainty of 3% in the overall tracking efficiency (±3%). 

• track reconstruction uncertainties close to the borders of acceptance This was 
evaluated by 

— raising (lowering) the cut on P^ mck to 160 MeV (140 MeV); 

— requiring \r)\ < 1.5 instead of 1.75; this effect dominates for Q 2 > 10000 GeV 2 ; 

— including tracks not associated to the primary vertex. 

• alignment uncertainties affecting the calculation of the boost vector to the Breit frame 
(^2%)- This was evaluated by 

— varying the polar angle for the scattered electron by ±2 mrad; 

— varying the polar angles for the hadrons by ±4 mrad. 

• assumptions concerning the details of the simulation of the hadronic final state (—4%). 
This was estimated by using Lepto instead of Ariadne. 

For data with Q 2 < 160 GeV 2 , the systematic uncertainties are slightly different and are 
due to: 

• imperfections in the simulation affecting the determination of the efficiency of event 
reconstruction and event selection 

• track reconstruction uncertainties (1q 5 %). 

• assumptions concerning the details of the simulation of the hadronic final state (+7%). 

• an uncertainty about the size of the contribution of diffractive events (1 2 %). This was 
estimated by varying the r) max cut by ±0.2 units. 

Further details can be found elsewhere [32]. All individual uncertainties were added in 
quadrature. 

6 Results 

6.1 Scaled momentum spectra 

Scaled momentum spectra were measured in the current region in the Breit frame as a 
function of Q 2 in the kinematic range 160 < Q 2 < 40960 GeV 2 and 0.002 < x < 0.75. 
Results are presented in Figs. HHH and Tables [2HHJ Also shown in Figs. dHU are previously 
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published results for 10 < Q 2 < 160 GeV 2 . In addition, in Fig. Q] data from a previous 
ZEUS publication [2] are given for 160 < Q 2 < 320 GeV 2 . They agree well with the 
measurements presented here. The same is true for previously obtained results up to 
Q 2 = 5120 GeV 2 which are not shown. 

The normalised spectrum, 1/N dn^/d ln(l/x p ), with N being the number of events and 
n* 1 being the number of charged particles, is shown in Figs. [IH2J These scaled momentum 
spectra exhibit a hump-backed form with an approximately Gaussian shape around the 
peak. The mean charged multiplicities are given by the integrals of the spectra. As Q 2 
increases, the multiplicity increases and, in addition, the peak of the spectrum moves to 
larger values of ln(l/x p ). 

In Fig. [U the predictions of Ariadne and Lepto are compared to the data. They 
reproduce the main features of the data but do not agree in detail. For the highest Q 2 
bin, both models predict too many charged particles at medium and low values of ln(l/x p ). 
Lepto also predicts too many particles for medium-Q 2 bins while Ariadne predicts too 
few for low-Q 2 bins. 

In Fig. [21 the MLLA+LPHD predictions [6], [12] are compared to the data. Too many 
particles are predicted for the highest- and lowest-Q 2 bins, while at medium Q 2 the data 
is reasonably well described. At low Q 2 , the observed particle deficit can be interpreted 
as a significant migration of particles to the target region of the Breit frame; this was 
also previously observed pQ. At medium Q 2 , the agreement is surprising, because BGF 
contributes significantly to the cross section and the observed particles should reflect the 
qq final state which is not included in the prediction. At the highest Q 2 available, the 
failure of the MLLA prediction probably reflects the fact that the e + e~ data used to 
obtain the input parameters are dominated by Z° exchange while, in ep collisions, photon 
exchange still dominates. 

The MLLA+LPHD calculations predict long tails towards large values of ln(l/x p ) over 
the complete range of Q 2 . These tails are sensitive to the mass correction applied in 
the calculation [Tallin]- The data do not show such tails in the predicted size. A better 
description of the large ln(l/x p ) region is obtained if m c g = Qo = 0.9A e g is taken instead 
of m cS = Q = A cff . 

6.2 Scaling violation 

As the energy scale, Q, increases, the phase space for soft gluon radiation increases, 
leading to a rise of the number of soft particles with small x p . These scaling violations 
can be seen when the data are plotted in bins of function of Q 2 . Figure [3] and 

Tables EHH] show that the number of charged particles increases with Q 2 at low x p and 
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decreases with Q 2 at high x p . Neither Lepto nor Ariadne provides a good description 
of this Q 2 dependence over the whole range of x p . 

Figure 0] shows the data together with four NLO+FF QCD predictions [SHTTj for x p > 0.1, 
where theoretical uncertainties are small and the predictions not too strongly affected by 
hadron-mass effects which are not included in the calculations [M]. The fragmentation 
functions (FF) used in all four calculations were extracted from e + e~ data. The four 
predictions are similar in shape and have similar uncertainties. The uncertainties are 
only illustrated for the calculation of Kretzer [5]. The NLO calculations also do not 
provide a good description of the data. Too many particles are predicted at small x p and 
too few at large x p . In general, the scaling violations predicted are not strong enough. 

Figure shows the same data as Fig. H] together with results from HI [3] and from e + e~ 
experiments [H]. For a proper comparison, the the particle momenta from e + e~ data 
were scaled to half of the centre-of-mass energy as discussed in the introduction and the 
scale was set to Q = 2 Ei, eaxn , where -Ebeam is the beam energy. In addition, corrections for 
the different treatment of K° and A decays were applied. The overall agreement between 
the different data sets supports fragmentation universality. The presentation of the data 
using a linear scale as presented in Fig. [6] does, however, show some significant differences 
between e + e~ and ep, in particular around the Z° mass at 0.02 < x p < 0.2 and at low Q 2 
at 0.1 < x p < 0.2. 

6.3 Limiting fragmentation 

The concept of limiting fragmentation [16] is based on the assumption that a Lorentz- 
contracted object passes through another object at rest, leaving behind an excited state 
with properties depending neither on the energy nor the identity of the passing object. 
This excited state fragments into particles in a restricted window of rapidity, called the 
limiting-fragmentation region. In this region, the limiting-fragmentation hypothesis pre- 
dicts that the density of charged particles per unit of rapidity depends only on W. 

Limiting fragmentation has been observed in a variety of hadronic collisions [501 15T] . 
including nucleus- nucleus interactions [T7] . It was observed that in the region of limiting 
fragmentation the particle density increases linearly with the rapidity before reaching a 
plateau. The slope of the increase did not show a W dependence, but the height of the 
plateau increased with W. These features are illustrated in Fig. [71 Bialas and Jezabek [52] 
proposed a statistical model to explain the missing W dependence of the slopes. In this 
model, soft particle production in hadronic collisions is described in terms of multiple 
gluon exchanges between partons of the colliding hadrons and by the subsequent radiation 
of hadronic clusters. 
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The application of the limiting-fragmentation hypothesis to e + e~ annihilations is not 
straight-forward. However, again a behaviour as illustrated in Fig. [7] was observed, only 
in this case the slopes increase with W [5T|l53]. 

In the case of ep collisions, the passing object is the proton while the virtual photon 
exchanged in the interaction is the object "at rest". It is assumed to be the excited 
hadron which fragments [53] . 

Figures [HH9] and Tables [10l4l~5l present the density of charged particles per unit of pseu- 
dorapidity, r/ Brcit , for 10 < Q 2 < 10240 GeV 2 in bins of Q 2 and W as listed in Table HJ 
A region of linear rise and the onset of a plateau are observed in all bins. This supports 
the applicability of the hypothesis of limiting fragmentation to the case of ep collisions. 
However, for low Q 2 , the plateau is only reached in the target region, ^ Breit > 0. 

The slopes in the region of linear rise do not depend significantly on either Q 2 or W, as 
also demonstrated in Figs. [10] and [TTJ The lack of a W dependence indicates that the 
model of Bialas and Jezabek is also applicable for ep collisions. 

Figures. [8] and also show predictions from Ariadne and Lepto. Overall, Ariadne 
provides reasonable predictions for the whole range in Q 2 and W. Lepto, however, 
predicts a sizeable increase in the height of the plateau with Q 2 and W which is not 
observed in the data. The predictions in the plateau region are sensitive to the input 
parameters used in the fragmentation functions. The usage of input parameters derived 
from SMC data [51] in Lepto results in the prediction of a softer spectrum, reflected in 
a charged-particle density of up to 30 % too high. 

The hypothesis of limiting fragmentation was further tested by studying the charged- 
particle densities in the rest frame of the fragmenting object, i.e. the virtual photon. The 
^Breit distributions were rebinned by shifting event by event all entries by lnQ/m^ , thus 
scaling the available energy to the pion mass. The resulting distributions are shown in 
Fig. [TTJ The distributions are very similar but for Q 2 > 160 GeV 2 a slightly larger slope 
is observed. This is a region where the BGF contribution is decreasing. In general, the 
observations support the hypothesis of limiting fragmentation and the model of Bialas 
and Jezabek. This indicates that, even at high Q 2 , soft processes are involved in the 
fragmentation and a statistical approach is justified. 

7 Conclusions 

Scaled momentum spectra have been measured in NC DIS for the current region in the 
Breit frame over the large range of Q 2 from 10 GeV 2 to 40960 GeV 2 . Large scaling viola- 
tions are observed. Comparing the data to e + e~ results generally supports the concept of 
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quark-fragmentation universality. Neither MLLA+LPHD nor NLO+FF calculations de- 
scribe the data well. A better, albeit not prefect description is provided by the Ariadne 
program. 

The limiting-fragmentation hypothesis and the statistical model of Bialas and Jezabek 
were tested by studying the density of charged particles as a function of the pseudorapidity, 
?7 Breit , over the range of 10 < Q 2 < 10240 GeV 2 . A region of linear rise and the onset of a 
plateau are observed over the whole range in Q 2 and support the limiting fragmentation 
hypothesis. The independence of the slopes on W supports the statistical approach of 
Bialas and Jezabek. 
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QIa (GeV 2 ) 


scaled momenta 


limiting fragmentation 




no. of events 


W (GeV) 


no. of events 


10 30 






80 160 


1 142 102 


30 160 






80 200 


1 185 843 


160 320 


0.0024 0.0500 


684 077 


80 240 


584 578 


320 640 
640 1280 


0.0100 0.0500 


170 784 
72 268 


120 240 
120 280 


152 367 
75 176 


1280 2560 


0.025 0.15 


31 608 


120 280 


29 949 


2560 5120 


0.05 0.25 


10 858 


160 280 


8 031 


5120 10240 
10240 20480 


0.05 0.50 


4 748 
1 197 


160 280 


2 914 


20480 40960 


0.05 0.75 


205 







Table 1: Number of accepted events in (Q ,x) bins for the scaled momentum analysis 
and in (W ,Q 2 ) bins for limiting-fragmentation studies. 
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]n(l/x p ) 


160 < Q 2 < 320 GeV 2 


320 < Q 2 < 640 GeV 2 


0.00 - 0.25 


0.0319 ± 0.0005j£ooi2 


0.0259 ± 0.0009±g;gg?g 


0.25 - 0.50 


0.1156 ± 0.001li°;°039 


0.101 ± 0.002+°;°°^ 


0.50 - 0.75 


0.2573 ± 0.0016±g;ggS 


0.232 ± 0.003±g:8o? 


0.75 - 1.00 


0.451 ± 0.002+°;°^ 


0.419 ± 0.004+°;°^ 


1.00 - 1.25 


0.695 ± 0.0031^21 


0.654 ± 0.005lg;g?| 


1.25 - 1.50 


0.959 ± 0.003+°;°™ 


0.925 ± 0.006+°;°^ 


1.50 - 1.75 


1.216 ± 0.004+°;°™ 


1.191 ± 0.007±g;g|| 


1.75 - 2.00 


1.437 ± 0.004+°;°^ 


1.440 ± 0.008ig;g^ 


2.00 - 2.25 


1.587 ± 0.004+°;°^ 


1.658 ± 0.008±g;g|| 


2.25 - 2.50 


1.643 ± 0.004+°;°|} 


1.772 ± 0.008ig;^f 


2.50 - 2.75 


1.599 ± 0.004+°;°™ 


1.861 ± 0.008+°;°| 


2.75 - 3.00 


1.459 ± 0.004+°;°^ 


1.829 ± 0.008±8:ob3 


3.00 - 3.25 


1.215 ± 0.003+°,;°,! 


1.673 ± 0.008^°,;°,™, 


3.25 - 3.50 


0.927 ± 0.003+^| 


1.398 ± 0.007+^ 


3.50 - 3.75 


0.649 ± 0.002+°;°^ 


1.087 ± 0.006^| 


3.75 - 4.00 


0.4209 ± 0.0019lg;S?7 


0.778 ± 0.005±g;g|l 


4.00 - 4.25 


0.2561 ± 0.0015+°;°°^ 


0.517 ± 0.004+°;°^ 


4.25 - 4.50 


0.1473 ± O.OOllj^oSti 


0.317 ± 0.003±g;gJ? 


4.50 - 4.75 


0.0842 ± 0.0009±S:88ii 


0.184 ± 0.002+°;°°j 


4.75 - 5.00 


0.0474 ± 0.0007ig;g^ 


0.1079 ± 0.0019iHoM 


5.00 - 5.25 


0.0266 ± 0.0006±°,;°,°,°J 


0.0598 ± 0.00151°,;°,°,^ 


5.25 - 5.50 


0.0167 ± 0.0006+°;°°°^ 


0.0356 ± 0.0013+°;°°^ 


5.50 - 5.75 


0.0102 ± 0.0006iHoi7 


0.0208 ± 0.0012+°;°°^ 


5.75 - 6.00 


0.0069 ± 0.0008^°,;°,°,°^ 


0.0126 ± 0.0012+°;°°^ 


6.00 - 6.25 


0.0040 ± o.ooni ,; , ,^ 


0.0092 ± 0.0016ig;gg^ 


6.25 - 6.50 


0.008 ± 0.007±g;gg| 


0.004 ± 0.002+°;°^ 



Table 2: The bin- averaged scaled momentum spectra, 1/N dn /dln(l/x p ), for 
160 < Q 2 < 640 GeV 2 . The first uncertainty is statistical, the second systematic. 
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ln(l/x p ) 


640 < Q 2 < 1280 GeV 2 


1280 < Q 2 < 2560 GeV 2 


0.00 - 0.25 


0.0228 ± 0.0013±£m 


0.0236 ± 0.0016+^3 


0.25 - 0.50 


0.086 ± 0.003l[!;^ 


0.087 ± 0.004+^;°^ 


0.50 - 0.75 


0.218 ± 0.004+Q oil 


0.199 ± 0.006+°oJe 


0.75 - 1.00 


0.371 ± 0.006+o q}i 


0.368 ± 0.009+00^ 


1.00 - 1.25 


0.597 ± 0.007+Q022 


0.577 ± 0.011+nmR 
" 1 -j- w.w j- j._o.016 


1.25 - 1.50 


0.857 ± 0.009+^^ 


0.840 ± 0.014^°-°^ 


1.50 - 1.75 


1.121 ± 0. 010+°,-°,? 7 

0.037 


1.098 ± 0.015+^^ 


1.75 - 2.00 


1.346 ± 0.011 + nnll 


1.325 ± 0.017+nn^Q 


2.00 - 2.25 


1.598 ± 0.012+n°Jn 
j-.u^w 0.050 


1.535 ± 0.017+nn^ 


2.25 - 2.50 


1.780 ± 0.013+nnco 
luu " " - 0.058 


1.779 ± 0.018±o°5? 


2.50 - 2.75 


1.920 ± 0.013+nn^ 
^ " " " - 0.058 


1.928 ± 0.019+nn^ 
j-.^i-^ 0.056 


2.75 - 3.00 


2.027 ± 0.013+2-2S 


2.135 ± 0.019^-^ 2 


3.00 - 3.25 


1.995 ± 0.013lg;^ 


2.152 ± 0.019+^11 


3.25 - 3.50 


1.816 ± 0.012 + n n^ 

0.055 


2.121 ± 0.018+nn^ 
0.063 


3.50 - 3.75 


1.575 ± 0.011 + nn«« 

0.068 


1.953 ± 0.017+nnl« 
w w 1 —0.058 


3.75 - 4.00 


1.241 ± 0.009+°,°i°, 

u-""^— 0.037 


1.705 ± 0.016+nn^ 
1 ww - 1 — 0.052 


4.00 - 4.25 


0.901 ± 0.008+^^ 


1.412 ± 0.014+°'°^ 

- 1 - ^ ^ —0.042 


4.25 - 4.50 


0.603 ± 0.006^^22 


1.038 ± 0.012 + S-So| 


4.50 - 4.75 


0.376 ± 0.005^;^ 


0.731 ± 0.0101^33 


4.75 - 5.00 


0.232 ± 0.004+°;°° 2 7 


0.463 ± 0.007+^;^ 


5.00 - 5.25 


0.133 ± 0.003±g;gg| 


0.288 ± 0.006±g;g?| 


5.25 - 5.50 


0.070 ± 0.002+°;°° 2 


0.163 ± 0.005±°;°°J 


5.50 - 5.75 


0.047 ± 0.002±g;gg| 


0.097 ± 0.004+°;°°^ 


5.75 - 6.00 


0.023 ± 0.002+°;°° 2 


0.049 ± 0.003±g;gg§ 


6.00 - 6.25 


0.013 ± 0.002+°;°° 2 ; 


0.026 ± 0.0031°,;°$ 


6.25 - 6.50 


0.007 ± 0.003±g;g?g 


0.017 ± 0.004+°;°°^ 


6.50 - 6.75 


0.024 ± 0.022+°;°°? 


0.011 ± 0.006+^ 



Table 3: The bin-averaged scaled momentum spectra, 1/N dn^ /d\n(l/x p ), for 
640 < Q 2 < 2560 GeV 2 . The first uncertainty is statistical, the second systematic. 
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ln(l/x p ) 


2560 < Q 2 < 5120 GeV 2 


5120 < Q 2 < 10240 GeV 2 


0.00 - 0.25 


0.0246 ± 0.0028±g:ggf? 


0.0261 ± 0.0044+°;°°^ 


0.25 - 0.50 


0.105 ± 0.008±Ho3 


0.075 ± 0.010+°;°°* 


0.50 - 0.75 


0.210 ± o.oi2±{J:om 


o.i78 ± o.oi7lg;85t 


0.75 - 1.00 


0.362 ± 0.016lo:{jii 


0.362 ± 0.026l8;8to 


1.00 - 1.25 


0.58 ± 0.02+°;^ 


0.630 ± 0.03518:81 


1.25 - 1.50 


0.78 ± 0.02+°;^ 


o.7oo ± o.o37i8:8?8 


1.50 - 1.75 


1.06 ± 0.03±g;g| 


1.02 ± 0.04j££ 


1.75 - 2.00 


1.32 ± 0.03l8;8! 


1.28 ± 0.05+°;83 


2.00 - 2.25 


1.49 ± 0.03±g:52 


1.54 ± 0.0518$ 


2.25 - 2.50 


1.73 ± 0.03±g:8§ 


1.87 ± 0.06^06 


2.50 - 2.75 


1.95 ± 0.03lg;JJ 


2.01 ± o.o5i8:o6 


2.75 - 3.00 


2.20 ± OMtom 


2.16 ± 0.05±8:oe 


3.00 - 3.25 


2.22 ± 0.03+°;55 


2.35 ± 0.05+8;o? 


3.25 - 3.50 


2.25 ± 0.03±g;gf 


2.40 ± 0.05l8:o? 


3.50 - 3.75 


2.27 ± 0.03^ 


2.53 ± 0.05l^| 


3.75 - 4.00 


2.13 ± 0.03lo:o7 


2.48 ± 0.05l^| 


4.00 - 4.25 


1.90 ± 0.03^ 


2.39 ± 0.051^1 


4.25 - 4.50 


1.51 ± 0.02±g5g 


2.06 ± 0.04+8;o6 


4.50 - 4.75 


1.17 ± 0.02lg;8g 


1.78 ± 0.04+8;8e 


4.75 - 5.00 


0.828 ± 0.017t^43 


1.45 ± 0.04i8:S 


5.00 - 5.25 


0.547 ± 0.0141K 


0.97 ± o.o3i8:8? 


5.25 - 5.50 


0.345 ± 0.011+°;°^ 


0.67 ± 0.0218:8s 


5.50 - 5.75 


0.188 ± 0.009l8:8°>8 


0.388 ± 0.018l8;8sJ 


5.75 - 6.00 


0.126 ± 0.008+°;°^ 


0.269 ± 0.017l8;82? 


6.00 - 6.25 


0.077 ± 0.008lg:oo5 


0.126 ± 0.015l8;822 


6.25 - 6.50 


0.033 ± 0.007l8:88g 


0.112 ± 0.018l8;8?3 


6.50 - 6.75 


0.027 ± O.OlOl^ 


0.058 ± 0.019lg;8n 



Table 4: The bin-averaged scaled momentum spectra, 1/N dn /d\n(l/x p ), for 
2560 < Q 2 < 10240 GeV 2 . The first uncertainty is statistical, the second system- 
atic. 
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in ^ i / j 


10240 < O 2 < 20480 GeV 2 


n n n k 
u.u — u.o 


n nfio j. n ni 1 +0.007 

U.UOZ It U.Ull^QQQJ 


n 1 n 
u.o — i.u 


n oon 4- n nos+O-039 

U.ZZU It U.UZO_Q QQg 


i.u — 1.0 


n c;7n 4- n n/i 7+0. 087 

U.OrfU it U.U4i_gQ 1 y 


1 t; on 
1.0 — z.U 


1.12 ± U.U/_ 004 


on o k 
z.u — Z.O 


1 Q7 4- n (17+0.19 

l.O/ it U.U ( _Q Q4 


o p; q n 
z.o — o.u 


1 4- n n8+0- 2 o 

Z.1Z it U.Uo_Q0g 


3.0 - 3.5 


2.60 ± 0.081q 10 


3.5 - 4.0 


2.56 ± 0.07±8;5| 


4.0 - 4.5 


2.60 ± 0.07±8;5| 


4.5 - 5.0 


2.16 ± 0.06±g;5| 


5.0 - 5.5 


1.36 ± 0.05io:§6 


5.5 - 6.0 


0.66 ± 0.03±g:S* 


6.0 - 6.5 


0.25 ± 0.031^2 


6.5 - 7.0 


0.04 ± 0.021^2 



Table 5: The bin-averaged scaled momentum spectra, 1/N dn /d ln(l/x p ), for 
10240 < Q 2 < 20480 GeV 2 . The first uncertainty is statistical, the second system- 
atic. 



ln(l/x p ) 


20480 < Q 2 < 40960 GeV 2 


0.0 - 1.0 


0.13 ± 0.04±g;g| 


1.0 - 2.0 


0.85 ± O.lOlaJjj 


2.0 - 3.0 


1.51 ± 0.12+^ 


3.0 - 4.0 


2.31 ± 0.13±8:S 


4.0 - 5.0 


2.43 ± 0.11±g;^ 


5.0 - 6.0 


1.32 ± 0.07±g;lS 


6.0 - 7.0 


0.33 ± 0.051^5 



Table 6: The bin-averaged scaled momentum spectra, 1/N dn* 1 /d Ya{l/x p ), for 
20480 < Q 2 < 40960 GeV 2 . The first uncertainty is statistical, the second system- 
atic. 
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< Q 2 >, GeV 2 


< x p < 0.02 


218 


9.32 ± 0.03±g;^ 


440 


i9.io d= o.o9±8:|| 


871 


35.53 ±0.181};$ 


1767 


6i.3 ± o.3i|:5 


3530 


93.6 ± 0.711$ 


6870 


140.1 ± 1.3lf;J 


13380 


191 ± 31J 1 


25700 


218 ± 7tu 


< Q 2 >, GeV 2 


0.02 < x p < 0.05 


218 


25.86 ± OMtoZ 


440 


39.34 ±0.111};^ 


871 


52.15 ±0.18^;| 


1767 


61.6 ±0.3l 2 J 


3530 


68.3 ± 0.5l|^ 


6870 


74.6 ± 0.8j£Jj 


13380 


78.6 ± 1.7l|;f 


25700 


69 ± 4+1 


< Q 2 >, GeV 2 


0.05 < x p < 0.1 


218 


19.79 ± OMtoH 


440 


25.34 ± 0.071^74 


871 


26.67 ±0.1llg;g 


1767 


27.20 ±0.1618:79 


3530 


27.6 ± 0.3lJ:l 


6870 


28.0 ±0.518:° 


13380 


26.8 ± 0.918:? 


25700 


23.3 ± 2.ll^:| 



Table 7: The number of charged particles per event and unit ofx p , 1/N n / Ax p , 
as a function ofQ 2 in bins of x p with widths Ax p . The first uncertainty is statistical, 
the second systematic. 
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< Q 2 >, GeV 2 


0.1 < x p < 0.2 


218 


9.338 ± 0.015ig;fj] 


440 


10.39 ± 0.03±g;|g 


871 


9.92 ± 0.05l^° 


1767 


9.68 ± 0.07lg;|| 


3530 


9.51±0.13l°;i 


6870 


9.6±0.2l°j 


13380 


8.5 ± OAtH 


25700 


7.9 ± 0.915J 


< Q 2 >, GeV 2 


0.2 < x p < 0.3 


218 


4.03 ± o.oiijj:^ 


440 


3.89 ± omtoit 


871 


3.60 ± omtoil 


1767 


3.52 ± 0.04l°i8 


3530 


3.29 ± 0.08±g:ig 


6870 


3.06 ±0.1218:^ 


13380 


2.87 ± 0.24±° o f 8 


25700 


3.1±0.6lgj 


< Q 2 >, GeV 2 


0.3 < x p < 0.4 


218 


1.806 ± 0.007l^g 


440 


1.717 ± 0.0131°;™ 


871 


1.515 ±0.01918;^ 


1767 


1.49 ± 0.03l°,o4 


3530 


1.51 ±0.05lHi 


6870 


1.65±0.09l8;o4 


13380 


1.33±0.16l8;i2 


25700 


o.7±o.3l8;4 



Table 8: The number of charged particles per event and unit of x p , 1/N n^/Axp, 
as a function ofQ 2 in bins ofx p with widths Ax p . The first uncertainty is statistical, 
the second systematic. 
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< Q 2 >, GeV 2 


0.4 < x p < 0.5 


218 


0.868 ± 0.0051S 


440 


0.785 ± 0.0091S 


871 


0.716 ±0.013±°;° 22 


1767 


0.721 ± 0.019±g;g|? 


3530 


0.692 ± 0.036±g;g|? 


6870 


0.714 ± 0.057±g;gf| 


13380 


0.392 ± 0.086±g;gf§ 


25700 


0.68 ± 0.34±g;i| 


< Q 2 >, GeV 2 


0.5 < x p < 0.7 


218 


0.329 ± 0.002±g;gig 


440 


0.295 ± 0.004±g;gg| 


871 


0.269 ± 0.006lg;g^ 


1767 


0.249 ± 0.008+°;°^ 


3530 


0.277 ±0.0151^08 


6870 


0.197 ± 0.02012:^6 


13380 


0.257 ±0.045lg;g^ 


25700 


0.33 ±0.15+^ 


< Q 2 >, GeV 2 


0.7 < x p < 1.0 


218 


0.056 ± 0.001±g;gg| 


440 


0.046 ± 0.001±g;ggi 


871 


0.041 ± 0.002+°;°^ 


1767 


0.042 ± 0.002±g;gg| 


3530 


0.049 ± 0.004±g;ggf 


6870 


0.043 ± 0.006±g;gg| 


13380 


0.054 ± 0.012±g;gg| 


25700 


0.013 ±0.013i^ 2 



Table 9: The number of charged particles per event and unit ofx p , 1/N n / Ax p , 
as a function ofQ 2 in bins of x p with widths Ax p . The first uncertainty is statistical, 
the second systematic. 
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^Breit 


160 <W< 200 GeV 


200 < W < 240 GeV 


240 < W < 280 GeV 


-5.5 -5.2 


0.15 ± 0.11 


0.054 ± 0.027 tofoe 


0.017 ± 0.012 toml 


-5.2 -4.9 


0.14 ± 0.08 t^t 


0.04 ± 0.02 t°o°ot 


0.027 ± 0.016 toml 


-4.9 -4.6 


0.07 ± 0.03 t°o°ot 


0.08 ± 0.03 t o°ot 


0.068 ± 0.024 toitl 


-4.6 -4.3 


0.15 ± 0.04 toil 


0.22 ± 0.04 tow 


0.16 ± 0.03 toil 


-4.3 -4.0 


0.25 ± 0.05 toil 


0.20 ± 0.04 tom 


0.23 ± 0.04 toil 


-4.0 -3.7 


0.33 ± 0.06 toil 


0.42 ± 0.06 t^ol 


0.28 ± 0.04 toil 


-3.7 -3.4 


0.51 ± 0.07 to.lt 


0.61 ± 0.06 t° 2 ot 


0.54 ± 0.05 tlit 


-3.4 -3.1 


0.97 ± 0.09 toil 


0.92 ± 0.08 toil 


0.64 ± 0.06 toil 


-3.1 -2.8 


1.31 ± o.io toil 


1.07 ± 0.08 toil 


0.96 ± 0.07 toil 


-2.8 -2.5 


1.54 ± 0.10 t°of 7 


1.37 ± 0.08 toil 


1.36 ± 0.08 toil 


-2.5 -2.2 


1.88 ± 0.10 tof 5 


1.60 ± 0.09 toil 


1.72 ± 0.09 ±g;g 


-2.2 -1.9 


2.05 ± 0.10 toll 


2.01 ± 0.09 toil 


2.09 ± 0.09 toil 


-1.9 -1.6 


2.34 ± 0.10 toil 


2.28 ± 0.09 toil 


2.46 ± 0.10 t°of 7 


-1.6 -1.3 


2.42 ± 0.10 toil 


2.47 ± 0.09 toil 


2.54 ± 0.10 toil 


-1.3 1.0 


2.30 ± 0.09 toil 


2.63 ± 0.09 toll 


2.81 ± 0.10 toil 


-1.0 -0.7 


2.52 ± 0.09 t°of 4 


2.82 ± 0.09 toil 


2.98 ± 0.11 t°;g 


-0.7 -0.4 


2.67 ± 0.09 toil 


2.68 ± 0.09 to.ll 


2.94 ± 0.11 tlil 


-0.4 -0.1 


2.41 ± 0.09 toil 


2.66 ± 0.09 to™ 


2.82 ± 0.10 toit 


-0.1 0.2 




2.46 ± 0.08 tom 


2.81 ± 0.10 tlil 


0.2 0.5 




2.81 ± 0.10 toil 


2.80 ± 0.10 tlil 


0.5 0.8 






3.08 ± 0.11 lg;| 



Table 10: The charged-particle density, 1/N dn± / drf reit , for 5120 < Q 2 < 
10240 GeV 2 and 5" bins in W . The first uncertainty is statistical, the second 
systematic. 
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^Breit 


160 < W < 200 GeV 


200 < W < 240 GeV 


240 < W < 280 GeV 


-4.6 -4.3 


0.112 ± 0.019 to°otl 


0.074 ± 0.015 ±°;°28 


0.069 ± 0.014 toiti 


-4.3 -4.0 


0.16 ± 0.02 t° m 


0.15 ± 0.02 tom 


0.095 ± 0.015 toml 


-4.0 -3.7 


0.24 ± 0.03 toil 


0.19 ± 0.02 ±g;g| 


0.19 ± 0.02 tof 3 


-3.7 -3.4 


0.35 ± 0.03 t° f 4 


0.34 ± 0.03 toil 


0.30 ± 0.03 toil 


-3.4 -3.1 


0.71 ± 0.04 toii 


0.57 ± 0.04 toil 


0.52 ± 0.03 toil 


-3.1 -2.8 


0.86 ± 0.04 tii) 


0.94 ± 0.04 toil 


0.70 ± 0.04 toil 


-2.8 -2.5 


1.31 ± 0.05 toll 


1.16 ± 0.05 tom 


1.08 ± 0.05 toil 


-2.5 -2.2 


1.59 ± 0.05 toil 


1.60 ± 0.05 tow 


1.43 ± 0.05 tUt 


-2.2 -1.9 


1.78 ± 0.05 t io 


1.80 ± 0.06 toil 


1.72 ± 0.05 toir 


-1.9 -1.6 


2.15 ± 0.06 toil 


2.11 ± 0.06 to^l 


2.11 ± 0.06 toil 


-1.6 -1.3 


2.28 ± 0.06 toil 


2.37 ± 0.06 toil 


2.32 ± 0.06 toil 


-1.3 -1.0 


2.28 ± 0.05 toil 


2.50 ± 0.06 toil 


2.45 ± 0.06 toir 


-1.0 -0.7 


2.30 ± 0.05 toil 


2.64 ± 0.06 toil 


2.63 ± 0.06 toi 2 7 


-0.7 -0.4 


2.46 ± 0.05 t°oi 7 


2.62 ± 0.06 toil 


2.85 ± 0.07 ±g;|| 


-0.4 -0.1 


2.41 ± 0.05 toll 


2.77 ± 0.06 toil 


2.82 ± 0.06 toil 


-0.1 0.2 




2.55 ± 0.06 toil 


2.66 ± 0.06 toil 


0.2 0.5 




2.77 ± 0.06 tof 5 


2.86 ± 0.06 toil 


0.5 0.8 




2.71 ± 0.07 toll 


3.16 ± 0.07 toil 


0.8 1.1 






2.92 ± 0.07 



Table 11: The charged-particle density, 1/N dn^ / dr] Breit , for 2560 < Q 2 < 
5120 GeV 2 and 3 bins in W . The first uncertainty is statistical, the second sys- 
tematic. 
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^Breit 
/ 


120 < W < 160 GeV 


160 < W < 200 GeV 


200 < W < 240 GeV 


240 < W < 280 GeV 


t: . O t: . W 


103 + 010 + ' 052 
w • i v 7 o _i_ w . w x w nni4- 


071 -1- 008 + 013 

u.u i x _i_ w . w wo — fl 01 


071 -1- 008 + ' 013 

W . W 1 X _l_ W . W WO a mo 


057 + 007 + ' 022 

W . WO 1 _l_ W.WWI n ni 1 


_4 n - -3 7 

t: . w o . I 


170 + 013 

u.i ( u zn u.u±o _o.036 


161 + 01 2 + 2 ™? 
u.iui _i_ u.uii -0.025 


127 + 010 

w . x zj i _i_ w . w x w 027 


088 + 008 

W.WOO _I_ W.WWO n ni n 


-3 7 - -3 4 

O.I O . t: 


2fi8 + 01 fi 

W.ZjWO —1— W « W X W n nQ^_ 


264 -1- 015 +™o? 

W « Zj Wt: _l_ U.UlU rt 0*34. 


218 + 013 

U.i 1U _l_ W.WXO n nofj 


203 + 012 

W.ZjWO —1— W.WXZj n nil 


_3 4 - -3 1 


53 + 02 + H2 

V 7 . '.J O _l_ W . W Zj n Qpr 


414 + 018 + ^SZ2 

u.<±_L4: _i_ u.uio -0.032 


350 + 01 7 
u.ooj in u.ui ( _o.035 


300 + 014 
u.ouu in u.u±<± _o.035 


-3 1 - -2 8 

O . X _ . o 


73 + 02 

W. 1 O _l_ W . W n no 


65 + 02 

u.UU _l_ W.WZj n n^ 


61 + 02 

U.UI _l_ U.Ui ri rj^ 


524 ±0 010 

W . OZjt: _1_ U.ui J n n9Q 


-2.8 - -2.5 


1 02 + 03 

1 .UZj _1_ W.WO n nc 


93 + 03 

W. <_/0 _l_ U.U<J n nc; 


81 ± 02 

\J . O X _l_ U.Ui n no 


77 + 02 

W. 1 1 _l_ W.WZj n 


-2.5 - -2.2 


1 41 + 03 +0 - 20 
±.^±i zc u.uo _ 009 


1 28 + 03 +0 - 05 

X . ZjO _1_ W . W O n ny 


1 10 + 03 +0 - 10 
i.lj zn u.uo _ 006 


1 00 + 03 +0 ' 07 

X . W 'J _l_ W.WO n nn 


-2.2 - -1.9 


1 67 + 03 +H 3 . 
±.u ( in u.uo _o.09 


1 62 + 03 

X.WZj _l_ W.Ul) n QjS 


1 50 + 03 + £ ™ 

X.WW _l_ W.WtJ n nn 


1 42 + 03 
i.-iz, zc u.uo „o.08 


-19 - -16 

x . '.y x . w 


1 93 + 03 

X.^yO _ 1— W.WO n nn 


1 02 + 03 

X . (7Zj _l_ W.UU n 19 


1 85 + 03 + £™ 

X . O O _l_ v ^ .WO n no 


1 74 ± 03 

X . 1 t: _l_ W.WO 10 


-1 fi - -13 

-L . W X . O 


2 09 + 03 + 249 


2 11 + 03 + S'S!o 
z,.±± in u.uo „ .08 


2 08 + 03 

Zj.WO —1— W.WO n 11 


2 04 + 03 + JH? 

Z.Ut: _l_ W.WO n 1 c 


-1 3 - -ID 

X . O J- . w 


2 21 + 03 + S'?o 

Zj.ZjX _l_ W . W O n io 


2 33 + 03 


2 26 + 04 

— — W _l_ W.W^ n ny 


2 25 + 04 

Zj.ZjO —1— W . W T: n 10 


-1 (1 - -0 7 

X . W W.I 


2 24 + 03 +0 - n 

Zj . Zjt: _l_ W.WO n nn 


2 42 + 03 +0 ' 10 

— • -L— _l_ U.UU n nn 


9 41 + 04 +0 - 10 

Zj . t: X _1_ W.W^ n 19 


2 45 ± 04 +0 ' 24 

Zj . tiO _l_ W . Wt: n 1^ 


-0 7 - -0 4 

W.I 1 7 . t: 


2 28 + 03 +242 

Zj.ZjU _l_ W.WO n in 


2 53 + 03 

^•Uu _l_ u.uiJ n 11 


2 60 + 04 + S4? 

Zj.WW —1— W . W t: n 10 


2 62 + 04 

Zj.WZj _l_ W.W^t n 10 


-0.4 -0.1 


2.31 ± 0.03 + ^-|? 

— U. -L 1 


2.65 ± 0.03 +0 r, -li 


2.64 ± 0.04 toil 


2.68 ± 0.04 toil 


-0.1 0.2 




2.60 ± 0.03 toll 


2.66 ± 0.04 tom 


2.68 ± 0.04 toil 


0.2 0.5 




2.54 ± 0.03 tofo 


2.81 ± 0.04 toil 


2.95 ± 0.04 toil 


0.5 0.8 




2.43 ± 0.04 toil 


2.52 ± 0.04 


3.08 ± 0.04 


0.8 1.1 






2.54 ± 0.04 toil 


2.88 ± 0.04 toil 


1.1 1.3 








2.79 ± 0.04 toil 



Table 12: The charged-particle density, 1/N dn* / dr] Bieit , for 1280 < Q 2 < 
2560 GeV 2 and 4 bins in W . The first uncertainty is statistical, the second sys- 
tematic. 



^Breit 


120 < W : 160 GeV 


160 < W < 200 GeV 


200 < W < 240 GeV 


240 < W < 280 GeV 


-4.0 -3.7 


0.086 ± 0.006 ±g:Sxo 


0.087 ± 0.006 tofu 


0.073 ± 0.005 i^os 


0.069 ± 0.005 toml 


-3.7 -3.4 


0.185 ± 0.008 t ° Q fs 


0.159 ± 0.007 t f Q7 


0.138 ± 0.007 t ol 7 


0.130 ± 0.007 t°om\ 


-3.4 -3.1 


0.301 ± 0.010 ±g:Sxs 


0.263 ± 0.009 t f 3 


0.238 ± 0.009 t° °om 


0.224 ± 0.009 tows 


-3.1 -2.8 


0.476 ± 0.013 to° u 


0.444 ± 0.012 ±g:g®| 


0.387 ± 0.011 t^26 


0.364 ± 0.012 toml 


-2.8 -2.5 


0.735 ± 0.015 t° ° f 2 


0.683 ± 0.015 tomt 


0.617 ± 0.014 tomr 


0.559 ± 0.014 toml 


-2.5 -2.2 


1.080 ± 0.018 t° fi 


0.950 ± 0.017 tomt 


0.881 ± 0.017 to° lt 


0.874 ± 0.018 ±o!o35 


-2.2 -1.9 


1.418 ± 0.019 tftjg 


1.283 ± 0.019 ±° ° il 


1.25 ± 0.02 +j»g 


1.22 ± 0.02 toll 


-1.9 -1.6 


1.74 ± 0.02 t° f 6 


1.62 ± 0.02 tof 7 


1.58 ± 0.02 ±g;5§ 


1.54 ± 0.02 t°oll 


-1.6 -1.3 


2.03 ± 0.02 toll 


1.92 ± 0.02 t°oll 


1.90 ± 0.02 ±g;i? 


1.86 ± 0.02 toll 


-1.3 -1.0 


2.17 ± 0.02 t° f 7 


2.15 ± 0.02 toll 


2.i3 ± o.o2 ig:5? 


2.15 ± 0.03 toll 


-1.0 -0.7 


2.27 ± 0.02 


2.29 ± 0.02 to°o 


2.33 ± 0.03 toll 


2.31 ± 0.03 toil 


-0.7 -0.4 


2.31 ± 0.02 toll 


2.43 ± 0.02 toll 


2.48 ± 0.03 toll 


2.43 ± 0.03 


-0.4 -0.1 


2.37 ± 0.02 tom 


2.48 ± 0.02 tom 


2.56 ± 0.03 toll 


2.50 ± 0.03 toil 


-0.1 0.2 


2.45 ± 0.02 t° f 7 


2.55 ± 0.02 


2.59 ± 0.03 jftn 


2.59 ± 0.03 toil 


0.2 0.5 


2.35 ± 0.02 t°of 8 


2.51 ± 0.02 toit, 


2.77 ± 0.03 toil 


2.78 ± 0.03 toil 


0.5 0.8 




2.49 ± 0.02 toil 


2.58 ± 0.02 toit 


2.73 ± 0.03 toil 


0.8 1.1 






2.57 ± 0.03 toil 


2.63 ± 0.03 Igiie 


1.1 1.4 






2.59 ± 0.03 toil 


2.64 ± 0.03 toil 



Table 13: The charged-particle density, l/N dn^ / d7] Bieit , for 640 < Q 2 < 
1280 GeV 2 and 4 bins in W . The first uncertainty is statistical, the second sys- 
tematic. 



^Breit 


120 < W < 160 GeV 


160 < W < 200 GeV 


200 < W < 240 GeV 


-3.7 -3.4 


0.094 ± 0.003 t ° n 


0.083 ± 0.003 ±° ™l 


0.082 ± 0.004 to™ 


-3.4 -3.1 


0.168 ± 0.005 tg;g$ 


0.145 ± 0.004 toml 


0.134 ± 0.004 to^l 


-3.1 -2.8 


0.285 ± 0.006 ±oIoiI 


0.246 ± 0.006 toml 


0.241 ± 0.006 to oli 


-2.8 -2.5 


0.451 ± 0.007 toml 


0.404 ± 0.007 to° \l 


0.395 ± 0.007 toml 


-2.5 -2.2 


0.714 ± 0.009 toml 


0.646 ± 0.009 t° ml 


0.621 ± 0.009 tomt 


-2.2 -1.9 


1.022 ± 0.010 t° il 


0.969 ± 0.011 to° ll 


0.915 ± 0.011 


-1.9 -1.6 


1.384 ± 0.012 t Zt 


1.335 ± 0.012 ±g;g|| 


1.253 ± 0.013 ±g:gff 


-1.6 -1.3 


1.702 ± 0.013 t° ;lw 


1.652 ± 0.013 i^77 


1.618 ± 0.014 ±°;J°| 


-1.3 -1.0 


1.998 ± 0.013 toafi 


1.936 ± 0.014 tofej 


1.903 ± 0.015 


-1.0 -0.7 


2.192 ± 0.013 t° .lf 


2.166 ± 0.015 ±%.om 


2.129 ± 0.016 tof 6 l 


-0.7 -0.4 


2.305 ± 0.013 t° ^l 


2.350 ± 0.015 to,mt 


2.323 ± 0.017 1^68 


-0.4 -0.1 


2.356 ± 0.013 t° f£ 


2.434 ± 0.016 to.mi 


2.422 ± 0.018 tomt 


-0.1 0.2 


2.447 ± 0.013 toml 


2.536 ± 0.015 t^.oil 


2.522 ± 0.017 toilt 


0.2 0.5 


2.364 ± 0.013 "±0:0?? 


2.549 ± 0.015 t° f 9 t 


2.680 ± 0.017 t°o 2 om 


0.5 0.8 


2.332 ± 0.014 ±° -l™ 


2.505 ± 0.015 ijjirajj 


2.617 ± 0.017 tofo°2 


0.8 1.1 




2.468 ± 0.015 tofsl 


2.591 ± 0.016 ±g:g|S 


1.1 1.4 




2.500 ± 0.016 tofgl 


2.569 ± 0.017 tofol 



Table 14: The charged-particle density, 1/N dn± / dr] Brcit , for 320 < Q 2 < 
640 GeV 2 and 3 bins in W . The first uncertainty is statistical, the second system- 
atic. 
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^Breit 


80 < W < 120 GeV 


120 < W < 160 GeV 


160 < W < 200 GeV 


200 < W = 240 GeV 


-3.4 -3.1 


0.107 ± 0.002 tofol 


0.092 ± 0.002 tofot 


0.082 ± 0.002 jftojjjj 


0.079 ± 0.002 1°;°05 


-3.1 -2.8 


0.186 ± 0.003 ±g:S2o 


0.155 ± 0.002 


0.144 ± 0.003 tH?2 


0.136 ± 0.003 ±g:g?i 


-2.8 -2.5 


0.304 ± 0.004 t ° n 


0.270 ± 0.003 ±g;gj§ 


0.248 ± 0.003 ±g;gj| 


0.233 ± 0.004 t ° w 


-2.5 -2.2 


0.512 ± 0.005 t° w 


0.436 ± 0.004 


0.411 ± 0.004 tofu 


0.386 ± 0.005 ±g;g|? 


-2.2 -1.9 


0.790 ± 0.006 t° f 3 3 


0.695 ± 0.005 t° -° o ll 


0.646 ± 0.005 ±g:gio 


0.609 ± 0.006 lg;g? 


-1.9 -1.6 


1.119 ± 0.006 t°o il 


0.992 ± 0.006 to° ll 


0.959 ± 0.006 tof 3 t 


0.902 ± o.oo7 ±g:g|i 


-1.6 -1.3 


1.479 ± 0.007 ±g:oss 


1.349 ± 0.007 t°o lo 


1.288 ± 0.007 to° tl 


1.242 ± 0.008 ±g;{£g 


-1.3 -1.0 


1.783 ± 0.007 tolU 


1.679 ± 0.007 toltl 


1.619 ± 0.008 toZl 


1.568 ± 0.009 ±g:gl? 


-1.0 -0.7 


1.989 ± 0.007 toml 


1.947 ± 0.008 1^052 


1.918 ± 0.009 t° f 5 i 


1.880 ± 0.010 ±^079 


-0.7 -0.4 


2.126 ± 0.007 tolll 


2.140 ± 0.008 tof 6 t 


2.134 ± 0.009 tolto 


2.120 ± 0.011 toml 


-0.4 -0.1 


2.233 ± 0.007 t° fi 


2.277 ± 0.008 toltl 


2.294 ± 0.010 tofi 


2.288 ± 0.011 ±oIoef 


-0.1 0.2 


2.309 ± 0.007 t° 2 f 6 


2.421 ± 0.008 tofel 


2.417 ± 0.010 t° ™ 


2.424 ± 0.011 ±g;oro 


0.2 0.5 


2.190 ± 0.008 toM7 


2.376 ± 0.008 toml 


2.513 ± 0.010 jftjjg 


2.617 ± 0.011 t° 2 ml 


0.5 0.8 




2.354 ± 0.008 ±g:ggf 


2.485 ± 0.009 ±g;5g| 


2.587 ± 0.011 1°;^ 


0.8 1.1 




2.377 ± 0.009 to.wl 


2.501 ± 0.009 ^056 


2.569 ± 0.011 


1.1 1.4 






2.500 ± 0.010 ±oIoto 


2.587 ± 0.011 ±g:g?i 



Table 15: The charged-particle density, l/N dri^ / 'dr) BTeit , for 160 < Q 2 < 
320 GeV 2 and 4 bins in W . The first uncertainty is statistical, the second system- 
atic. 



^Breit 


80 < W : 120 GcV 


120 < W < 160 GeV 


160 < W < 200 GeV 


-4.0 -3.7 


0.034 ± 0.002+°;^ 


0.033 ± 0.002+°;°°^ 


0.027 ± 0.002+°;°°? 


-3.7 -3.4 


0.057 ± 0.002lg;gg? 


0.053 ± 0.003+°;°°^ 


0.049 ± 0.003+°;°°^ 


-3.4 -3.1 


0.093 ± 0.003j£ooi 


0.086 ± 0.003±g;gg? 


0.084 ± 0.0041°;°°^ 


-3.1 -2.8 


0.158 ± 0.004+°;°°! 


0.141 ± 0.0041°;°°/; 


0.138 ± 0.004+°;°°^ 


-2.8 -2.5 


0.254 ± 0.005±g;gg| 


0.248 ± 0.005^;^ 


0.214 ± 0.005±g;ggt 


-2.5 -2.2 


0.413 ± 0.006+^;°,^ 


0.361 ± 0.006±g;gJi 


0.339 ± 0.006+^;°^ 


-2.2 -1.9 


0.600 ± 0.007+°;°^ 


0.542 ± 0.007+°;°^ 


0.530 ± 0.008+°;°°^ 


-1.9 -1.6 


0.856 ± 0.0081^ 


0.796 ± 0.009±g;g^ 


0.773 ± 0.009j#8i2 


-1.6 -1.3 


1.140 ± 0.009+^ 


1.084 ± O.OlOl&jjn 


1.034 ± O.OlllHn 


-1.3 -1.0 


1.424 ± O.OlOlg;^ 


1.383 ± 0.011±g;gl| 


1.324 ± 0.012+°;°?o 


-1.0 -0.7 


1.682 ± O.Olllg;^ 


1.634 ± 0.012+°;°^ 


1.600 ± 0.013±g;gf| 


-0.7 -0.4 


1.855 ± O.Olllgim? 


1.833 ± 0.012+°;°^ 


1.827 ± 0.014+°;°^ 


-0.4 -0.1 


1.979 ± 0.012+°;°^ 


1.986 ± 0.013+°;°| 


2.039 ± 0.014+°;°^ 


-0.1 0.2 


2.060 ± 0.0121°;°^ 


2.134 ± 0.013j#o22 


2.142 ± 0.014±°,;°JJ 


0.2 0.5 


2.051 ± 0.012+°;°^ 


2.201 ± 0.013±g;^? 


2.301 ± 0.014+°;°| 


0.5 0.8 


2.082 ± 0.012+°+°° 


2.263 ± 0.013+^ 


2.327 ± 0.014+°+^ 


0.8 1.1 




2.294 ± 0.0131^ 


2.395 ± 0.0141°+^ 


1.1 1.4 




2.302 ± 0.013+°+°^ 


2.440 ± 0.015+°+! 


1.4 1.7 




2.287 ± 0.0141°;°^ 


2.497 ± 0.015±g:Sf| 



Table 16: The charged-particle density, 1/N dn ± / drj Bveit , for 60 < Q 2 < 
160 GeV 2 and 3 bins in W . The first uncertainty is statistical, the second system- 
atic. 
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Brcit 


80 < W < 120 GeV 


120 < W < 160 GeV 


160 < W < 200 GeV 


-3.1 


-2.8 


0.070 ± 0.005^;^ 


0.051 ± 0.004+°;°°^ 


0.057 ± 0.005±g:gg§ 


-2.8 


-2.5 


0.105 ± 0.005lg;^ 


0.096 ± 0.006^;^ 


0.073 ± 0.005±g:gg| 


-2.5 


-2.2 


0.177 ± 0. 007+0:007 


0.157 ± 0.007+o;oog 


0.139 ± O.008l8;8o7 


-2.2 


-1.9 


0.298 ± 0.009+°;°}^ 


0.271 ± 0.010lg;gio 


0.250 ± O.OlOlffi 


-1.9 


-1.6 


0.469 ± 0.0111°;^ 


0.436 ± 0.012±g;gl§ 


0.391 ± 0.0121°;°^ 


-1.6 


-1.3 


0.706 ± 0.0141°;°^ 


0.654 ± 0.0141°;°^ 


0.611 ± 0.016lH?g 


-1.3 


-1.0 


0.999 ± 0.016^;^ 


0.951 ± 0.017lg;gJ 


0.809 ± 0.017^.^6 


-1.0 


-0.7 


1.293 ± 0.018ig;g^ 


1.233 ± 0.019l{j;{jP 


1.17 ± 0.02+°;°^ 


-0.7 


-0.4 


1.534 ± 0.019^;^ 


1.51 ± 0.02±g;™ 


1.44 ± 0.02lg;{g 


-0.4 


-0.1 


1.80 ± 0.021°;^ 


1.75 ± 0.021°;^ 


1.68 ± 0.03+°;°^ 


-0.1 


0.2 


1.97 ± 0.02±g;g| 


1.96 ± 0.02±g;g| 


1.93 ± 0.03t{j:o3 


0.2 


0.5 


2.06 ± 0.02±g;g| 


2.07 ± 0.0218;™ 


2.19 ± 0.03±g;g| 


0.5 


0.8 


2.11 ± 0.02lg;J2 


2.17 ± omt^l 


2.24 ± 0.031^ 


0.8 


1.1 


2.17 ± 0.02±g;S2 


2.28 ± 0.0218:51 


2.30 ± omtoil 


1.1 


1.4 


2.24 ± 0.03±g;g| 


2.33 ± Omtoll 


2.48 ± 0.03l{J;$ 


1.4 


1.7 


2.18 ± 0.03±g;5§ 


2.36 ± 0.03lg;Jg 


2.56 ± 0.03lg;S| 



Table 17: The charged-particle density, 1/N dn ± / dr] Breit , for 50 < Q 2 < 60 GeV 2 
and 3 bins in W . The first uncertainty is statistical, the second systematic. 
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Brcit 


80 < W < 120 GeV 


120 < W < 160 GeV 


160 < W < 200 GeV 


-3.1 


-2.8 


0.049 ± 0.003+^;°,°^ 


0.044 ± 0.003±g;gg| 


0.039 ± 0.003+^;°^ 


-2.8 


-2.5 


0.083 ± 0.004±°;°°| 


0.089 ± 0.005±g;g^ 


0.064 ± 0.004^°;°°* 


-2.5 


-2.2 


0.149 ± 0. OOS^ooe 


0.136 ± 0.006+g;oo5 


0.118 ± 0.006+o;oog 


-2.2 


-1.9 


0.242 ± 0.006^°;°°* 


0.242 ± 0.008+°;°°j 


0.194 ± 0.007+°;°^ 


-1.9 


-1.6 


0.405 ± 0.0081°;°}^ 


0.384 ± 0.009+°;°^ 


o.33o ± o.oo9±g:g5| 


-1.6 


-1.3 


0.607 ± 0.010ig;g^ 


0.576 ± 0.011±g;g?l 


0.495 ± O.OllinS 


-1.3 


-1.0 


0.912 ± 0.013+^ 


0.815 ± 0.013±g;g?| 


0.766 ± 0.014iH?4 


-1.0 


-0.7 


1.176 ± 0.014±°;°^ 


1.121 ± 0.015^;°22 


1.029 ± 0.016±g;g?i 


-0.7 


-0.4 


1.461 ± 0.015+^ 


1.416 ± 0.017j£ra7 


1.39 ± 0.02+°;°^ 


-0.4 


-0.1 


1.710 ± 0.017+°;°™ 


1.649 ± 0.018±g;gff 


1.67 ± 0.02+°;°? 


-0.1 


0.2 


1.901 ± 0.018±g;8|^ 


1.883 ± 0.019ig;§26 


1.88 ± 0.021°;°^ 


0.2 


0.5 


2.009 ± 0.018±g;g|| 


2.044 ± 0.019±g;g| 


2.08 ± 0.021°;°^ 


0.5 


0.8 


2.106 ± 0.018±g:g|f 


2.14 ± 0.02+°;°* 


2.25 ± 0.02±°+° 


0.8 


1.1 


2.187 ± 0.018±g;g^ 


2.26 ± 0.02+°;^ 


2.30 ± 0.02±°;^ 


1.1 


1.4 


2.24 ± 0.02+°;°^ 


2.33 ± 0.02+°;^ 


2.39 ± 0.02+°+° 


1.4 


1.7 


2.21 ± 0.02±g;g| 


2.40 ± 0.02±g;o4 


2.46 ± 0.02±g;S| 



Table 18: The charged-parUcle density, 1/N dn^ / dr) Bleit , for 40 < Q 2 < 50 GeV 2 
and 3 bins in W . The first uncertainty is statistical, the second systematic. 
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^ Brcit 


80 < W < 120 GeV 


120 < W < 160 GeV 


160 < W < 200 GeV 


-3.1 -2.8 


0.036 ± 0.002i{j:ora 


0.035 ± 0.002l°:882 


0.029 ± 0.00218:883 


-2.8 -2.5 


0.066 ± 0.003i{j:oo| 


0.055 ± 0.00318:882 


0.056 ± 0.00318:884 


-2.5 -2.2 


0.122 ± 0.004+o;°o4 


0.108 ± o.oo4i8:88I 


0.090 ± 0.00418:883 


-2.2 -1.9 


0.196 ± 0.005±g;gg| 


o.i85 ± 0.00518:88? 


0.160 ± 0.005l8:88e 


-1.9 -1.6 


0.344 ± o.oo6i8;88| 


0.300 ± 0.00618:884 


0.260 ± 0.006l8:8oe 


-1.6 -1.3 


0.529 ± 0.007l8:8os 


0.488 ± 0.00818:81? 


0.431 ± 0.00818:889 


-1.3 -1.0 


0.767 ± 0.009l8:8?o 


0.696 ± 0.00918:808 


o.653 ± o.oio!8:8?8 


-1.0 -0.7 


1.050 ± o.oio!8;8u 


0.964 ± 0.011±8;g5* 


0.913 ± 0.01218:816 


-0.7 -0.4 


1.342 ± o.oni8:8o8 


1.260 ± o.oi2l8:8i? 


1.228 ± 0.01418:8m 


-0.4 -0.1 


1.597 ± 0.01218:8?? 


1.548 ± 0.01418:8^ 


1.508 ± o.oi6i8:8?2 


-0.1 0.2 


i.8i6 ± o.oi3i8:8I 


1.795 ± 0.01418:81 


1.757 ± o.oi6i8:8?I 


0.2 0.5 


1.973 ± 0.014l8:8ig 


i.98o ± o.oi5i8:8?4 


2.035 ± 0.01718:8™ 


0.5 0.8 


2.089 ± 0.01418:824 


2.133 ± 0.01518:826 


2.168 ± 0.01718:826 


0.8 1.1 


2.200 ± o.oi4i8:81 


2.229 ± 0.015l8:o?8 


2.303 ± 0.018l8:o32 


1.1 1.4 


2.268 ± 0.01518:8a? 


2.328 ± 0.016l8:o?2 


2.367 ± 0.018l8:o5? 


1.4 1.7 


2.285 ± 0.01718:8s? 


2.379 ± 0.016l8:o| 


2.473 ± 0.019l8:o32 



Table 19: The charged-particle density, 1/N dn ± / dr] Breit , for 30 < Q 2 < 40 GeV 2 
and 3 bins in W . The first uncertainty is statistical, the second systematic. 
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^Brcit 


80 < W < 120 GeV 


120 < W < 160 GeV 


-3.1 -2.8 


0.0268 ± 0.0012^;°°^ 


0.0233 ± 0.0012±g;ggi| 


-2.8 -2.5 


0.0475 ± 0.001618;^ 


0.0470 ± 0.0019±g;g^g 


-2.5 - -2.2 


0.091 ± 0.002+o;°o3 


0.075 ± 0.002+°;oo2 


-2.2 -1.9 


0.155 ± 0.003±g:8SI 


0.130 ± 0.003+°;°^ 


-1.9 -1.6 


0.256 ± 0.004±g:ooi 


0.220 ± 0.004+°;°^ 


-1.6 -1.3 


0.403 ± 0.0051^09 


0.364 ± 0.005±8:8o° 


-1.3 -1.0 


0.620 ± 0.006ig;gg? 


0.543 ± 0.006lg;gJ| 


-1.0 -0.7 


0.863 ± 0.007±g;g?| 


o.783 ± o.oo7l8;gg? 


-0.7 -0.4 


1.142 ± 0.008ig;g^ 


1.055 ± 0.008^07 


-0.4 -0.1 


1.427 ± 0.009+°;°}^ 


1.370 ± 0.009l8;8o7 


-0.1 0.2 


1.692 ± 0.009iHis 


1.634 ± 0.010l8;[!o 7 


0.2 0.5 


1.882 ± O.OlOlg;^ 


i.87i ± 0.01018:^2 


0.5 0.8 


2.049 ± 0.010±g:gJ| 


2.072 ± o.oni8:8I 7 


0.8 1.1 


2.177 ± o.oio±g:g|i 


2.265 ± 0.012l8:o23 


1.1 1.4 


2.269 ± 0.011±g;g^ 


2.365 ± 0.012+8:ol5 


1.4 1.7 


2.290 ± 0.011±g;g?| 


2.429 ± 0.012l8:oi 



Table 20: The charged-particle density, 1 /N dn* / dri BltAt , for 20 < Q 2 < 30 GeV 2 
and 2 bins in W . The first uncertainty is statistical, the second systematic. 
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^Brcit 


80 < W < 120 GeV 


120 < W < 160 GeV 


-3.1 


-2.8 


0.0156 ± 0.0006 ±olooo? 


0.0138 ± 0.0006^;^ 


-2.8 


-2.5 


0.0300 ± 0.0008^;^ 


0.0252 ± 0.0008^;^ 


-2.5 


-2.2 


0.0548 ± 0.001 1^:0020 


0.0458 ± 0.001 lloioow 


-2.2 


-1.9 


0.0915 ± 0.0014+°;°°^ 


0.0785 ± 0.0014±J™ 


-1.9 


-1.6 


0.1564 ± 0.0019+°;°°^ 


0.1327 ± 0.0018^;^4 


-1.6 


-1.3 


0.260 ± 0.002±g:gg| 


0.223 ± 0.002+°;°^ 


-1.3 


-1.0 


0.404 ± 0.003ig;gg3 


0.354 ± 0.003±g;88| 


-1.0 


-0.7 


0.591 ± 0.003^;^ 


0.529 ± 0.004+°;°^ 


-0.7 


-0.4 


0.827 ± 0.004±g;g?| 


0.769 ± 0.004ig;g^ 


-0.4 


-0.1 


1.100 ± 0.005±g:8ii 


1.051 ± 0.005±°;°^ 


-0.1 


0.2 


1.409 ± 0.005iHie 


1.369 ± 0.006±g;g|| 


0.2 


0.5 


1.665 ± 0.006±g;g^ 


1.667 ± 0.006^;^o 


0.5 


0.8 


1.900 ± 0.006±g;gfg 


1.929 ± 0.007lg;^ 


0.8 


1.1 


2.086 ± 0.007lg;g^ 


2.132 ± 0.0071^?? 


1.1 


1.4 


2.202 ± 0.007±g:oi5 


2.304 ± 0.008lg;J?J 


1.4 


1.7 


2.259 ± 0.007±g;g|i 


2.395 ± 0.008lg;S^ 



Table 21: The charged-particle density, 1/N dn ± / dr] Breit , for 10 < Q 2 < 20 GeV 2 
and 2 bins in W . The first uncertainty is statistical, the second systematic. 
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Figure 1: The scaled momentum spectra, 1/N dn^ /d\n{l/ x p ) , for different 
(x, Q 2 ) bins. The dots represent the new, the triangles the previous ZEUS mea- 
surement. The data overlap in the 160 < Q 2 < 320 GeV 2 bin. The inner error 
bars, where visible, indicate statistical uncertainties, the outer statistical and sys- 
tematic uncertainties added in quadrature. The full and dashed lines represent the 
Lepto and the Ariadne predictions, respectively. 
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Figure 2: The scaled momentum spectra, 1/N dn ± /d ln(l/x p ), for different 
(x,Q 2 ) bins. The band represents the range of the MLLA+LPHD predictions. 
Other details as in Fig. [IJ 
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Figure 3: The number of charged particles per event per unit of x p , 1/N n / Ax p , 
as a function of Q 2 in x p bins of width Ax p . Other details as in Fig. [IJ 
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Figure 4: TTie number of charged particles per event per unit of x p , 1/N n / 'Ax p , 
as a function of Q 2 in x p bins with width Ax p as in Fig. The shaded band repre- 
sents the NLO calculation by Kretzer /2J/ with its renormalisation scale uncertainty. 
Additional NLO calculations are shown: Kniehl, Kramer, Potter f^(KKP), Albino, 
Kniehl, Kramer JJ^(AKK) and De Florian, Sassot and Stratmann ITIf(DSS). 
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Figure 5: The number of charged particles per event per unit of x p , 1/N n / Ax p , 
as a function of Q in x p bins with width Ax p . Also shown are data from HI J3J/ 
and e + e~ [L$ . The dots (triangles) represent the new (previous) ZEUS measure- 
ment, the squares the HI data and the inverted triangles the e + e~ data. The inner 
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systematic uncertainties added in quadrature. The three lowest x p bins are scaled 
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Figure 6: The number of charged particles per event per unit of x p , 1/N n /Ax p , 
as a function of Q 2 in x p bins with width Ax p . Other details as in Fig. [5j 
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Figure 7: The main features of the prediction of the Bialas-Jezabek model based 
on the limiting-fragmentation hypothesis for the dependence of the particles density 
on the rapidity in hadronic collisions for two values ofW with Wl > W2. 
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Figure 8: The normalised charged-particle density per unit of ry Breit ; 

1/N dn ± /d?7 Breit ; for different {W,Q 2 ) bins for Q 2 > 160 GeV 2 . The dots rep- 
resent the ZEUS measurement. The error bars, where visible, indicate statistical 
and systematic uncertainties added in quadrature. The dashed and dotted lines 
represent the Lepto and Ariadne predictions, respectively. 
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Figure 9: 

1/N c/n ± /rfr/ Breit ; 
details as in Fig. 
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Figure 10: The normalised charged- particle density per unit of rj at , 

1/N c/n ± /(i?7 Breit ; for 4 bins in Q 2 and 5 bins in W . The error bars, where visible, 
indicate statistical and systematic uncertainties added in quadrature. 
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Figure 11: The normalised charged-particle density per unit of r] Brcit , 

1/N dn^ I 'dr] Bmit , for 2 bins in W and 5 bins in Q 2 . The distributions were re- 
binned by shifting the entries for each event by Y=lnQ/m 7T as explained in the 
text. The error bars, where visible, indicate statistical and systematic uncertainties 
added in quadrature. 
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